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Abstract 

We address in a recent gauge model of unparticles the issues that are important for consis- 
tency of a gauge theory, i.e., unitarity and Ward identity of physical amplitudes. We find that 
non-integrable singularities arise in physical quantities like cross section and decay rate from 
gauge interactions of unparticles. We also show that Ward identity is violated due to the lack of 
a dispersion relation for charged unparticles although the Ward-Takahashi identity for general 
Green functions is incorporated in the model. A previous observation that the unparticle's (with 
scaling dimension d) contribution to the gauge boson self-energy is a factor (2 — d) of the par- 
ticle's has been extended to the Green function of triple gauge bosons. This (2 — d) rule may 
be generally true for any point Green functions of gauge bosons. This implies that the model 
would be trivial even as one that mimics certain dynamical effects on gauge bosons in which 
unparticles serve as an interpolating field. 
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1 Introduction 



As the era of the Large Hadron Collider is approaching, many new theoretical ideas have been 
contemplated that could potentially be tested there. One of radical suggestions is perhaps that 
of unparticle by Georgi |[T]|. Such an object is by definition not a particle, but some stuff that 
follows scale invariance, though it may well arise from certain high energy scale theory of par- 
ticles. The scale invariance makes a dispersion relation generally not possible for an unparticle; 
instead, it determines its kinematics in terms of a parameter, called scaling dimension, of its 
corresponding field. The very nature of the invariance also implies that the field is generically 
non-local. The latter results in novel features not seen in the particle world, for instance, non- 
trivial interference in the time-like regime t2j (see also Ref t3J), one particle to one unparticle 
transitions [4|, and non-integral power laws of long distance forces between particles mediated 
by unparticles [5] (see also Ref [6J), etc. 

Unparticles must interact with particles to be physically relevant since we manipulate par- 
ticles in experiments, and the interactions can be systematically organized in effective field 
theory. Although most studies, both phenomenological and theoretical, cope with bosonic un- 
particles that couple as a standard model singlet to particles [|7]|, it is hard to imagine that un- 
particles must not carry the standard model charges. As a matter of fact, the first gauge model 
of unparticles has been constructed in Ref [[8]|. In this circumstance, fermionic unparticles 
ll9l[T0l[IIl[T2l[T3l can equally well couple to particles, and their phenomenology could be even 
more interesting than their bosonic counterparts lfT2ll . 

In this work we continue our theoretical investigation on the gauge model of |[8l and address 
some issues that have only been lightly touched upon in [14J. A gauge model of unparticles must 
pass the standard consistency criteria like unitarity and Ward identities for scattering amplitudes. 
We make these checks and find the answer is negative for both. This means that the model is 
not yet amenable to computing physical amplitudes involving unparticles in the initial or final 
state. In [[T4|. we also observed that the scalar unparticle contribution to the complete (not just 
the imaginary part of as shown in [8|) gauge boson self-energy is exactly (2 — d) times that of 
a scalar particle in the same representation, where d is the scaling dimension of the unparticle. 
We extend this to the case of triple gauge bosons. This seems to indicate that this (2 — d) rule 
is generally true. In that case, the model of |[8l would be naive even as one that mimics certain 
dynamical effects on gauge bosons in which unparticles serve as an interpolating field. 

The paper is organized as follows. We describe briefly the gauge model in the next section 
and catalog the Feynman rules to be employed in later sections. The derivation of the rules is 
outlined in the appendix. In section [3] we show explicitly that the (2 — d) rule holds true for 
the Green function of triple gauge bosons. This is then followed in section |4] by comparative 
unitarity checks for ungauged and gauged unparticles using the simplest two-point Green func- 
tions of particle fields. Although the Green functions in the gauge model fulfil Ward-Takahashi 
identities, we demonstrate in section [5] that physical amplitudes involving unparticles and phys- 
ical gauge bosons do not satisfy Ward identities. We conclude with some remarks in the final 
section. 
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2 Gauge model and Feynman rules 



The scale symmetry of a scalar unparticle field of scaling dimension d demands its inverse 
propagator to be proportional to {—p^ — ze)^^"^, with p being its momentum [|2l|3l|. This is a 
non-integral power for a general real number d>l, and thus corresponds to a non-local field. 
The non-locality makes the conventional minimal gauging not work. Fortunately, a similar 
non-local problem was successfully dealt with some years ago in the context of reproducing 
low energy Goldstone dynamics from dynamical quarks in QCD lfT5l[T6l . The lesson has been 
recently applied to gauging unparticle fields in [jSl. 

A scalar unparticle multiplet may be coupled to gauge fields via the Wilson line. The 
action is (HI 

S = j d^xdS^\x)E{x-y)F{x,y), (1) 



F{x,y) = Pexp 



-igT" J^'l, dw^ 



^(y), (2) 



where P denotes path-ordering that effects on the group generators T" in the unparticle represen- 
tation, and g is the gauge coupling. i^^E{z) is the Fourier transform of the inverse propagator: 



E(z) = I -0ji'-""=B(p), (3) 



with A{d) a J-dependent constant not essential for our purpose. Note that putting an infrared 
cut-off in the propagator does not modify our subsequent conclusions. 

The action ([T]) contains gauge interactions that are quadratic in unparticle fields but involve 
gauge fields to an arbitrary order. There is no obstacle to derive their Feynman rules though the 
procedure rapidly becomes more and more involved as the number of gauge fields increases. 
Some details of it are given in the appendix. We list below the vertices up to three gauge fields 
that will be required in later sections. 

The Feynman rules for the vertices up to two gauge fields are known in [[8]|. The A^'^^^ 
vertex is (with ig to be attached on both sides) 

r"^{-p-q,p;q) = T%2p + q)^Ei [p-q), (5) 

where momenta before the semicolon are for unparticles and those after it for gauge bosons, 
with all momenta being incoming. The A^A^'^ vertex is (with ig^ to be attached on both 
sides and differing by a factor of i from |[8l) 

^fvi-p-qn:P;quq2) = giiv{T''J^}Ei{p;qn) 

+ T"T''{2p + q2)v{2p + q2 + qi2)^E2{p;qi2,q2) 

+ T''T"{2p + qi)^{2p + qi+qi2)vE2{p;qi2,qi), (6) 
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where qi2 = + etc, and the following notations are used, 

E{a + b) -E{a) 



Ei{a;b) 
E2{a;bi,b2) 



[a + bf-a^ ' 
Ei{a\bi) -Ei{a;b2) 



{a + biy~{a + b2) 
E2{a;bi,c) -E2{a;b2,c 
{a + bi)^-{a + b2) 

The notations are slightly improved over those in lfT4l to better display symmetry. 
Finally, the A^^A^pAy^ vertex is (with ig^ to be attached on both sides) 

^apr'^-P - ^123 ,p;qi,q2,q3) 
= T^'iT", T^jgap {qm + qn + '^p)yE2{p;qm,q\2) 
+ {r\T^]rg^p{q^^2p)yE2{p\qn3.q3) 

+ T''T^T''{qi23 + q23 + ^p)a{q23+q3 + ^p)p{q3 + ^p)rE3{p;qm,q23;q3) 

+ T''T''T''{qi23+q3l+^P)l5iq3l +q3+2p)aiq3 + 2p)YE3ip;qi23,q3i,q3) 

+ 2 perms. (8) 



3 {2 — d) rule for triple gauge bosons 

When the action in ([T]) is exponentiated and integrated in the path integral over the unparticle 
fields, we obtain an effective action of the gauge fields. It is in this sense that the gauge model 
discussed here is parallel to the non-local chiral quark model in low energy hadronic physics [[TSl 
[T6l . In the latter case, integration over chiral quarks with a momentum-dependent dynamical 
mass results in an effective action for the Goldstone bosons and external sources or gauge fields. 
The chiral quark fields serve as an interpolating field that mimics the strong dynamics of QCD 
as manifested in the low energy constants in chiral Lagrangian. It therefore sounds reasonable 
to expect that the unparticle fields in the gauge model should at least play a similar role in the 
context of certain new strong dynamics at a high energy scale. 

1 1 
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Figure 1. Diagrams contributing to triple gauge boson function. 

In the previous work [[T4|. we observed that the unparticle contribution to the gauge boson 
self-energy follows a simple rule; namely, it is (2 — d) times the contribution from scalar par- 
ticles in the same representation. Its possible impact on the running and unification of gauge 
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couplings was also studied. In this section we examine whether the rule applies to the Green 
function of triple gauge bosons which is kinematically more varied than a self-energy. If it does 
apply, it would unlikely be accidental but might be generally true. 

The contributing diagrams to the function are shown in Fig. 1, in which the double dashed 
line stands for the scalar unparticle and the three gauge bosons carry the group (a,Z7,c) and 
Lorentz indices (a, /3, 7) with the incoming momenta qi. The vertices in section [21 yield 

/d'^p 
j2^^ai-p + ^2, p + ^y^1l)^i-p - p;q3)'i^p{-p, p - q2;q2) 

X D{p)D{p + q3)D{p-q2) + {q2:li,b) ^ (^3,7,^), 

Kpf = -ig^irj -^r"^{-p-qi,p;qi)rf^{-p,p + qi-,q2,q3)D{p)D{p + qi) 

+ 2 perms, 

<pf = ishr j ^T%{-p,p;q,,q2,q,)D{p). (9) 

The integrals can be defined in n dimensions for regularization, but our subsequent algebraic 
manipulation does not depend on it. We stress again that introducing an infrared cut-off to the 
propagator does not modify our discussion either. The particle case can be recovered in the 
limit J — > 1 whence E\ ^ I: The graph (3) vanishes identically and the graph (2) vanishes due 
to symmetry, while the two terms in graph (1) combine to 



^3rabc^(. f A {2p + q3-q2)a{2p~q2)p{2p + q3)y 
particle ^ ^ j (2;r)4 [p2 + /g] [(p + ^73)2 + /g] [(p - ^j)^ + /c] ' ^ ^ 



where trT"!^ = C{r)5"^ for particles in representation r. 

In the unparticle case, the graph (3) contains two classes of terms. The integrand in the first 
(class I) is proportional to 

E,{p;Q) = limE,{p;q) = ~ (11) 

q^O p^ + l£ D[p) 

while the remaining terms constitute the class II whose integrand is not proportional to E\ {p;0). 
We demonstrate that the class II terms in graph (3) cancel completely the graphs (1) and (2). 
First, the terms in class II that involve a single signature tensor cancel in pair upon taking the 
traces and doing integration, so do the similar terms in graph (2). Second, by inspection, the 
remaining terms in II and graph (2) can be combined in pair. Choosing judiciously the routing 
momenta and making use of identities of Ei, we found that those terms condense to a sum of 
two 'form factors' with the common Lorentz structure, {2p — q2 + q-i)"'{2p + qi,y{2p — q2)^ • 
This is the same Lorentz structure exactly for the first term in -^^^^y^ and up to a minus sign 
for the second upon flipping the sign of p. The first form factor is (with the prefactors ig^ 
understood). 



dj 



dkei[ei-ej)^ 1 

rij-rik Ui - Uj 



(12) 
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where the sum is over the set (123), (231), (312). We have introduced some abbreviations: 



n\=p^, n2 = {p-q2f, n^, = {p + q2,Y ; 
d\=D{p), d2 = D{p-q2), dj, = D{p + q-i); 

(13) 

Uj-n^ 

where ijk is again cyclic in 123. The second form factor is obtained from the first by the 
interchanges of indices, b ^ c and 2^3, but with a global minus sign. Third, using identities 
of fractions, it is straightforward to show that the sum in eq (fT2l) is equal to Ji J2^3^i^2^3- This is 
again the integrand in the second term of -^^^y upon extracting the Lorentz structure displayed 
above so that the sum in (fT2l) (i.e., the first form factor) completely cancels the second term in 
"^(x^y ' Similar cancellation occurs between the second form factor and the first term in eq (fT2)) . 
To summarize, the complete graphs (1) and (2) are cancelled by class II terms in graph (3). 

Now we are left with class I terms in graph (3). First of all, eqs (II 1I9|) imply that the 
integrand is of a particle type. The terms with a signature tensor are again cancelled in pair. 
Using the fraction identity, ^ [ni{ni — nj) (n,- — H/t)]^ = [^i^2^3]~\ some algebra verifies our 

m 

claim: 



^abc 
^ajiy 



unparticle 



jabc 
^aPy 



(14) 

particle 



We end this section with a remark. It seems unlikely that the above relation is specific to 
two- and three-point functions of gauge fields. Our explicit demonstration of it might suggest a 
way to reach the general result for any point functions: The cancellation mechanism witnessed 
in two- and three-point functions might indicate that the only contribution for any point function 
comes exclusively from the tad-pole like graph involving the highest point vertex in each case. 



4 Unitarity 

Unitarity of the scattering matrix is one of the fundamental criteria that any quantum theory 
must meet. This is especially true of a gauge theory in which additional delicacies may occur. 
The purpose of this section is to show using the simplest possible process that the gauge model 
of unparticles proposed in [8J breaks unitarity. For comparison, we also examine unitarity in 
non-gauge interactions between unparticles and particles, and we find that these interactions 
generally preserve unitarity in the conventional sense. 

Consider the one-loop self-energy of a scalar particle field arising from interactions with 
scalar unparticles. To the graphs (1) and (2) shown in Fig. 2 there correspond the two effective 
interactions: 

ifi = Ai$(p^, (15) 
^2 = h^Wi, (16) 
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where <I>, (p are the scalar particle fields of mass M, m, and ^ , '^i, ^2 the scalar unparticle 
fields of scaling dimension J, Ji, d2, respectively. 



k 



„'-"" k " X- 



(1) (2) 
Figure 2. Self-energy of scalar particle field arising from eqs (I15I16I) . 

The imaginary part of the self-energy in graph (2) is found to be, for M > m. 



= (1 - rY'-' 2F,{d-Wd; 1 - r^)Bid,d), (17) 

where r = m/M and 2^1 and 5 are the standard special functions. This should be compared to 
the decay width for 4» — > 9 -|- for unitarity check. Finishing all phase space integrals but that 
of the unparticle energy yields 



■ (18) 



Changing the variable to m = — 5(1 +'')]/[5(l — works out the integral to 



X 2Fi|-^,J-l;J + yi:p^ 



5 J-1, 



(19) 



The unitarity relation Im £/i (M^) = MFi is verified using the relation 

F(^rf-l,J;2rf;^^)=22(^-i)(l+r)-2(^-i)F(-l J-l;J + l;z^^ 



(20) 



and a relation for B function. 

The interaction ^2 involves two unparticles and has been less discussed in the literature. Its 
contribution to the imaginary part of graph (2) is easily worked out to be 

Im^2(M2) = _ A|A(Ji)A(J2) ^2(..+.,-2) 
^ ^ 4(4;r)2 



X .I^^-^^:^B{dud2)- ^i"(^i+^2);r 



r(2-Ji)r(2-j2) 
7 



sin(Ji;r) sin((i2^) 



(21) 



The phase space for the decay 4> — '^i'^2 is more involved. Finishing integrals of one unparti- 
cle momentum and the angles of the other, we obtain 



where 



ll^dvodvvHvl-vY[{l-vo)- 



.21/? 



(22) 



(23) 



Here the integration region R in the vvq plane is bounded by the lines, v = 0, v = vq, and 
v + vo = 1. As the integrand is even in v, we make the region symmetric under v —v. The 
integrals are then factorized by the new variables, vq — v = x, vq + v = y with x G [0, 1] and 
y G [0, 1] so that, fora>— 1,/3>— 1 (i.e., Ji 2 > 1 in our case). 



2-' [' dx dy {y-x)^{xyn{l-y){\-x)f 
Jo Jo 



= 2 



-3 



5(a + 3,/3 + l)5(a + l,/3 + l)-(5(a + 2,/3 + l))^ 



r(a + l)r(a + 2)r(i8 + l)r(j8 +2) 



(24) 



8r(a + /3+3)r(a + /3+4) ' 

The unitarity relation = Im ■e^(M^) is confirmed using r(z)r(l — z) sin(z;r) = n. 

It is not surprising that unitarity is preserved by non-gauge interactions of unparticles be- 
cause the unparticle propagator has the correct cut structure by construction ^ and because 
those interactions are Hermitian. As we pointed out in [[T4||. the gauge interactions of unparti- 
cles in [m are actually non-Hermitian in the time-like regime. This may be a source of unitarity 
violation in the model. In Ref [[T4|. we reached this conclusion by symmetry analysis for the 
process, qq ^ '^'^ via gluon exchange assuming is charged under QCD. In what follows, 
we demonstrate the violation analytically by the simplest possible process of the gauge boson 
decay, A'^{p) ^ ^{ki)^ {k2), with the gauge boson momentum p in the time-like regime. 

The imaginary part of the gauge boson self-energy from the unparticle loop can easily be 
obtained from that of the scalar particle's using the (2 — d) rule: 



Im n(/) 



f PfiPv 



-g^LV ^{P 



(2-rf)4^C(.)p 



(25) 



The amplitude for the decay is basically the vertex shown in eq ©. The properly summed and 
averaged decay rate is 



g^C{r)^/p^sm^{d7t)J{d), 



where 



J{d) 



dvodv 



(l-2vo)^ 



(l-Vo)2-V^ 



2-d 



+ 



(1-vo)^ 



2-d 



(26) 



(27) 
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with the same region R as in r2 above. The terms in the square brackets arise from combination 
of phase space factors and the numerator of Ei in the vertex. The potential singularity at vq = 5 
is due to the denominator in Ei, and makes it impossible to factorize the integral into two 
separate ones as we did in 1^ . 

To observe the non-integrability that the singularity may cause, we finish the v integral first. 
Using the symmetry with respect to vq = 5, we restrict ourselves to the left half of R and obtain 

where 

f{a;t) - f'B(^,\+a^ 2Fi(^a,^;^+a;t 

+ t-"B(^^A-a^ 2Fi(^-a,^;^--a-,t^-^. (29) 

Note that J{d) is even in {2 — d) (so is F) while Im Tl{p^) is odd. This is the argument employed 
in [fT4| to signify unitarity breakdown. But what really occurs is even worse: The singularity 
introduced by the vertex in eq Q is logarithmically non-integrable. To see this, one needs to 
expand f{a; 1 — z) at z = 0: 

/(a;l-z) = a[\{f{\+a)-\{f{\-a)]z 

+ 0a^lnz + a^[C-v/(l-a)-VA(l+a)] 

7 7 \ 

--a{l-a)\j/{2-a) + -a{a +l)\j/{2 + a) j + 0{z^), (30) 

where = r'{^)/r{^) and C is a constant. The unitarity is thus badly violated in this case 
by non-integrable singularities introduced in interaction vertices of the gauge model. Note in 
passing that there is no problem with the particle limit of J ^ 1 although it is better to take the 
limit at the very start to avoid the ambiguity between the sin^( J;r) factor in T and the singularity. 
We stress that this breakdown of unitarity occurs at all energy scales in the gauge model of 
unparticles, in contrast to the conventional effective field theory in which unitarity starts to be 
violated at energy scales close to its ultraviolet cut-off. This also implies that unitarity cannot 
be simply recovered by including new degrees of freedom in the gauge model as we do in a 
conventional effective theory. 



5 Ward identity 

Now we address the issue of Ward identity necessary for a consistent gauge theory. We do not 
expect any problem with Ward-Takahashi identity for generally off- shell Green functions as it is 
built in by construction of the gauge model. As examples, we list below the first few identities 
for Green functions involving up to three gauge bosons. 
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The simplest one is 

q^Il{-p-q,p;q) = r' [D \p + q) - D-\p)], 
while the next one requires a little rearrangement of terms, 

(Ii^fvi-P-Qn,p;qhq2) 



(31) 



-p-gi2,p + qi;g2)T" -T'T''y{-p-q2,p;q2 



+ ir^'T'^{~p-qx2.P\qn). 



(32) 



The derivation of identity for the triple gauge boson vertex is much more involved. The main 
trick is to use partial fraction. But when the cloud clears up, the answer is simple: 



fobd Y^dc 



-p-qm,p;q\,q2,q3] 



{-p-qi23,p+qi;q2,q3)T" - TT'p^i-p - ^23 , ^2,^3 



+ ir'^'^{-p-qn3.P\ql2.q3)-ir''T%{-p-qm.P\q2.q^ 



(33) 



In the conventional gauge theory of particles, the Ward identity for physical amplitudes is 
obtained from the Ward-Takahashi identity by going to the physical limit of charged particles. 
In an Abelian theory like QED it is sufficient to require electrons to be on-shell. But in a non- 
Abelian theory like QCD, it is necessary that gluons be physical as well since they are also 
charged. In a gauge theory of unparticles however, an on-shell condition (dispersion relation) 
is missing for unparticles; this may endanger the Ward identity for physical amplitudes. If this 
happens, unphysical states of gauge bosons can be produced by unparticles, which is of course 
not acceptable. We show below by a simple process that this happens indeed in the considered 
model. 

Consider the process of unparticle pair production by the fusion of a gauge boson pair, 
A^(/:i)A^(fc2) — * ^ (Pi)'^ {P2), whose Feynman diagrams are depicted in Fig. 3. Putting the 
gauge bosons on-shell, kf =k2 = 0, the once contracted component amplitudes are 



ab;l 



X 



-[r , r^] [k,p{pl - pD + ^2/3^1 ■{P\-P2) + {P2 - Pl)l52ki ■ k2\ 

Ei{-pi;pi+P2) 



ab:2 



2k\ ■ k2 



{r,r}fci^£i(-pi;pi+p2) 

2pi +k2)p [Ei{-pi;pi+p2) -Ei{-pi;k2)] 



+ T^T''^^^{-2pi+2ki+k2)p [Ei{-pi;pi+p2) -Ei{-pvM)] ■ 



k2-p2 



ab;3 



T^T"{2p2-k2)p 



h -pi 



ab-A 



J-'Orpb 



k2-P2 

D{p2-ki) 



1- 



D{p2-k2) 



D{P2) 



D{P2) 

{-2p\ +k2)pEi{-pi;k2) 



Ei{-piM] 



(34) 
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where a factor is implied on the right hand side. The particle case is recovered by sending 
D^^{q) and E\ — > 1. As the above result seems hopeless, we examine the simpler, less 
restrictive, doubly contracted amplitude. The sum is 



bTCl 



+ 



D{pi)D{p2 



(^Dik2-p2)-^[D{pi)+Dip2)]^ , (35) 



which does not vanish as the Ward identity requires. The situation does not improve in the 
Abelian case. 





'\/\/\/\/\yr- 




(2) 



(3) 



(4) 



Figure 3. Diagrams contributing to the process A^(/:i)A^(/:2) — > ^{pi)'^iP2)- 

A simple way out seems to require D^^{p) = for unparticles appearing in the initial and 
final states. While the meaning of this is obscure by itself, it implies a dispersion relation for 
unparticles with d < 2. But this is obviously not a consistent concept as there would be no 
difference between a physical particle and unparticle. Furthermore, it would break unitarity 
established for unparticles that have no gauge interactions. 



6 Conclusion 

It looks natural that unparticles are charged under the standard model gauge group. But it is 
rather difficult to couple them to gauge fields since they are generically nonlocal in nature. 
Nevertheless, the first gauge model has been attempted in Ref [8] (for recent discussions about 
the work and development, see [[TtI [T8l [T9l ). We have investigated in this work whether the 
model fulfils the standard requirements that a consistent gauge theory must do; namely, the 
unitarity and the Ward identity for physical amplitudes involving unparticles. We find that the 
answers to both are negative. 

In non-gauge interactions no surprise is expected for unitarity as the unparticle propagator 
incorporates correct analyticity properties and the interactions are usually trivial in analyticity 
structure. The latter is no longer the case in the gauge model considered here. We find that 
its gauge interactions introduce non-integrable singularities in phase space that make physical 
quantities like cross section meaningless. This failure in unitarity occurs at any energy scale and 
thus cannot be cured by incorporating new degrees of freedom as one does in a conventional 
effective field theory. The result on Ward identity may be surprising at first sight since the 
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Ward-Takahashi identity for general Green functions has been built in the model. The point 
here is that to obtain the Ward identity for physical amplitudes some physical conditions have 
to be imposed to delete the contact terms in the Ward-Takahashi identity. These are the on-shell 
condition for charged particles and the transversality condition for gauge bosons if charged. It 
is the lack of a dispersion relation for charged unparticles that the passage is blocked. These two 
defects might easily be blamed upon the wisdom of conformal field theory that no consistent 
scattering matrix is known. But this does not explain why we seem to obtain reasonable results 
with non-gauge interactions of unparticles by following the standard procedure in field theory. 

Even if a gauge model of unparticles is afflicted with these diseases, it could still serve as a 
useful tool to mimic certain strong dynamical effects on gauge bosons at low energies. In that 
case, unparticles appear as an interpolating field confined to the virtual loops of gauge bosons. 
This situation is parallel to the relationship between the nonlocal chiral quark model lfT5l[T6l and 
the low energy dynamics of Goldstone bosons. Because of this, we have considered the Green 
function of triple gauge bosons due to unparticle loops and found that an earlier observation 
on the self-energy of gauge bosons also applies here. Namely, the unparticle contribution to 
the Green function is a factor (2 — d) of the scalar particle's in the same representation of the 
gauge group. We conjecture that this (2 — d) rule may be generally true. If this were the case, 
the model would be too naive even if the unparticles are considered as an interpolating field. It 
looks fair to say that the challenge of gauging unparticles still remains. 

Acknowledgement This work is supported in part by the grants NCET-06-0211 and NSFC- 
10775074. 

Appendix. Derivation of Feynman rules 

We outline the derivation of the vertices in eqs (I5I6I8I) for completeness. The basic technique 
was developed in Ref ilSl . The vertices with up to two gauge bosons were known previously 
ifTSl lSl for which our manipulation of series is slightly simpler, while the vertex with three gauge 
bosons is not yet available in the literature. Introducing Fourier transforms of the functions 




(36) 



and assuming E{p) is a function of (which is the case here), expansion dX p^ =0 yields for 
the action in eq ©: 





;i=0 
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Replacing e'P'^{p^)" = [{—d^)"e'P'y] and doing integration by parts, the action becomes 



n=o n. J I \ Jx J 

which is the starting point to all vertices. We have used the following abbreviations: 



(dk) 



{2ny 



dx = d X, d{x) = 5 (jc), 



and the derivatives always refer to y unless otherwise stated. 

We begin with the derivation of eq ([5]). The vertex in coordinate space is 

5^S 



5A«(xi)5^(v)5^t(2) 



- ldy5{y-z){-dY[i-is)LT5{y-v)] 

n=0 



where from now on the following short-cuts will be used, 

laa^ja 8{X,-U)du'', 5, = 5 {Xi - y) . 

Its Fourier transform yields 
£W(0) 



n=0 



dxidydze'^P'^'Py-'i'''^ {-d^y5{z-y) [-/^Lf] • 



The basic trick here is integration by parts. The n = term vanishes, while the n = I term j 
using d^L"^ = T'dig^", 

{n=\ term) = E^^\Q)gT''{p + pY {2nf5{p'-p-q). 
The general term is obtained by induction. Assuming 



(n— th term) 



-gr\p + prfn{2Kf8{p'-p-q), 



with /o = 0, /i = 1, the coefficient for the next term is found to be 

fn+l=P^fn + {p + qf\ 

with k^" = [k^Y; namely. 



/n+l 



{p+q) 



2n 



1+r 



fn 



(p+^)2("-l) 



1 - r"+^ 

l+r + --- + r" = — , r 

1 — r 



{p+qy 
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The general coefficient is thus 



fn 



Using 



£W(0) (p + g)2«_p2« ^ E{p + q)-E{p) 
' n\ [p + qY-p^ {p + q)^-p^ 



(47) 



(48) 



eq ([5]) is obtained. 

Now we derive the vertex in eq Q: 



J dxidx2dvdz 



5^5 



^(")(0) 



n=0 



dxidx2dydz \{-d^) "5{y- z) 



J- all J bv 
^1 ^2 



(49) 



The n = term again vanishes, while the n = I term gives 

(n=l)term = g2£(i)(0)gMV|7^«^7^^} (^27T)'d{p' - p - qn), 

with qi2 = ^1 + q2- Since not all possible Lorentz structures have appeared, we have to go one 
step further. The n = 2 term yields after some algebra 



{n = 2) term = ^ ' 



{lnf5{p' -p- qn) {g^'{T\ T^W^ + p' 



2! 

+T^T%2p + q,f{2p + q,+qnY 
+T^T'^{2p + q2y{2p + q2 + qi2r}. 



(50) 



Assuming the coefficients for the {T", T*}, T^T" and T^T^ terms to he gn{p,qi,q2), ^1,^2) 
and /?„(/?, ^2, <?i) respectively, with the initial conditions: 



(51) 



go = 0, gi = 1, g2 = p'^+P^; ho = hi =0, /Z2 = 1, 
we find out their (n + l)-th expressions after some work: 

gn+iip^Quqi) = p^gn{p,quq2) + {p+qi2f", 

hn+l[p,q2,q\) = P hn{p,q2,qi) + -, — ■ ^ — 7 — ■ r^-- 

[p + qnY - {p + q2Y 

gn has the same structure as /„ while /z„, when multiplied by [{p-\- qnY ~ + I2Y] , becomes 
a difference of two series each having the same structure as /„ again. We thus find 

{p + quY^'-p''" 



(52) 



gn{p,quq2) 

K{p,q2,qi) 



(p+qnY-p^ ' 
1 

{p+qi2Y-ip+q2Y 



[p + quY^-P^" ip + q2Y"-p^" 



{p+q 
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{p+q2) 



(53) 
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Then eq © obtains readily. 

Finally, we describe briefly the derivation of the vertex with three gauge bosons. The algebra 
blows up rapidly as the number of gauge bosons increases. The vertex to compute is 

5^S 



jndx.dvdze dAUxM 



. ^ J ^^^^^^^^ ,'[p'z-j:,,r^]s^y_^^ (_^2y^ |e-'P-^'P[LrLf Lf]} . (54) 



Both n = and n = 1 terms vanish. For n = 2 term, we use 



e-ipy [{pY-4p^ppd^dP +4ip^ppdP-2p^d^ 
-lip^d^d^-lippd^dP + (^2)2] p[Lf Lf Lf ], 



bfi J cy-i 



(55) 



and note that only terms with three or more derivatives can contribute because of the d{y — z) 
and that for the same reason only those without after differentiation survive. Extracting out 
{-i)ig^^E^'^\Q){2KY5{p' -p-l^qi), the result is 

r-{r,r*}(^i23+^i2+277)^ + {r,r^}r(^3 + 2p)^]g«^+2perms. (56) 

Since the non-^ terms have not appeared, we have to compute explicitly the n = 3 term. This is 
the most tedious part for the vertex as it involves derivatives up to the sixth order. One should 
be very careful that derivatives do not commute because of the path ordering operation. We 
skip the detail of the calculation but recording the result. Leaving aside the common factors 
(-0?'g^j!^^^HO)(2;r)45(y -7? -E^/), the n = 3 term is 

'(T^{T\T'}[p^ + {p + qnf + {p + qn,)\qx2, + qn + 2py 

r*}r-[/ + + ^3)2 + + ^j23)'] (^3 + 2pf) g"^ + 2 perms 
+ [rr^r(^i23 + ^23 + 2p)"(^23+^3 + 2;?)^te + 277)^ + 5perms]. (57) 

Now we assume the above structure is valid for the n-th term and denote the coefficients of the 
displayed three terms as bn{p,qi,q2,q3), Cnip^qi^qi^qi), dnip.qi^qi.qs)- The initial condi- 
tions are therefore 

b2 = C2=d3 = 1, 

^3(77,^1,^2,^3) = P^ + ip + qnf + ip + qmf, 

^3(^,^1,^2,^3) =P^ + iP + q3f + iP + qm)^: (58) 



while those not listed vanish. The coefficients in the next term are found to be (with arguments 
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p, ^1,^27 ^3 suppressed), 



{p+qmr-{p+qi2r 

2 , {p + qm)^"-{p + q3?" 



C„_|_i — p Cn + 



{p+qmY-ip + q^Y 



dn+l = p^dn + 



{p + qmY-{p + q23Y 
{p + ^123)'" - (P + q3Y" {P + qisY" -ip + qsY" ^ 

{p + qmY-{p + q3Y {p + qisY - {p + qsY 



(59) 



bn and c„ have the same structure as while d„, when multiplied by [{p + qnsY ~ (P + qisY], 
is a difference of two series each having the structure of /i„. They are worked out to yield the 
final answer shown in eq ([8]). 
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